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In this work we consider Randall- Sundrum brane- world type scenarios, in which the spacetime 
is described by a five - dimensional manifold with matter fields confined in a domain wall or three- 
brane. We present the results of a systematic analysis, using dynamical systems techniques, of 
the qualitative behaviour of Friedmann - Lemaitre - Robertson - Walker type models, whose matter 
is described by a scalar field with an exponential potential. We construct the state spaces for these 
models and discuss how their structure changes with respect to the general - relativistic case, in 
particular, what new critical points appear and their nature and the occurrence of bifurcation. 
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I. INTRODUCTION 

String and membrane theories are considered to be 
promising candidates for a unified theory of all forces 
and particles in Nature. Since a consistent construction 
of a quantum string theory is only possible in more than 
four spacetime dimensions we are compelled to compact- 
ify any extra spatial dimensions to a finite size or, alter- 
natively, find a mechanism to localize matter fields and 
gravity in a lower dimensional submanifold. 

Motivated by orbifold compactification of higher di- 
mensional string theories in particular the dimensional 
reduction of eleven - dimensional supergravity introduced 
by Hofava and Witten ||, [^], Randall and Sundrum 
showed that for non - factorisable geometries in five di- 
mensions there exists a single massless bound state con- 
fined in a domain wall or three - brane || . This bound 
state is the zero mode of the Kaluza- Klein dimensional 
reduction and corresponds to the four - dimensional gravi- 
ton. This scenario may be described by a model consist- 
ing of a five - dimensional Anti - de Sitter space (known 
as the bulk) with an embedded three - brane on which 
matter fields are confined and Newtonian gravity is effec- 
tively reproduced on large - scales. A discussion of earlier 
work on Kaluza - Klein dimensional reduction and matter 
localization in a four - dimensional manifold of a higher - 
dimensional non - compact spacetime can be found in |Q . 

Gravity on the brane can be described by the standard 
Einstein equations modified by two additional terms, one 
quadratic in the energy - momentum tensor and the other 
representing the electric part of the five - dimensional 
Weyl tensor. 

Following an approach to the study of homogeneous 
cosmological models with a cosmological constant first 
introduced by by Goliath and Ellis |B| (see also || 
for a detailed discussion of the dynamical systems ap- 
proach to cosmology), Campus and Sopuerta have re- 
cently studied the complete dynamics of Friedmann - 
Lemaitre - Robertson - Walker (FLRW) and the Bianchi 
type I and V cosmological models with a barotropic equa- 
tion of state, taking into account effects produced by the 
corrections to the Einstein Field Equations . 

Their analysis led to the discovery of new critical points 



corresponding to the Binetruy - Deffayet - Langlois (BDL) 
models 0, representing the dynamics at very high en- 
ergies, where effects due to the extra dimension become 
dominant. These solutions appear to be a generic feature 
of the state space of more general cosmological models. 
They also showed that the state space contains new bifur- 
cations, demonstrating how the dynamical character of 
some of the critical points changes relative to the general - 
relativistic case. Finally, they showed that for models 
satisfying all the ordinary energy conditions and causal- 
ity requirements the anisotropy is negligible near the ini- 
tial singularity, a result first demonstrated by Maartens 
et. al. §]. 

In this paper we extend their work to the case where 
the matter is described by a dynamical scalar field 4> with 
an exponential potential V(<t>) — exp(60). This work 
builds on earlier results obtained by Burd and Barrow 
|]lOf and Haliwell |n[] who considered the dynamics of 
these models in standard general relativity (GR). Van 
den Hoogcn et. al. and Copeland et. al. have also 
considered the GR dynamics of exponential potentials 
in the context of Scaling Solutions in FLRW spacetimes 
containing an additional barotropic perfect fluid |l2| , p"3[ . 

It is worth mentioning that Exponential potentials are 
somewhat more interesting in the brane -world scenario, 
since the high - energy corrections to the Friedmann equa- 
tion allow for inflation to take place with potentials or- 
dinarily too steep to sustain inflation E3. 



II. PRELIMINARIES 

In what follows we summarise the geometric framework 
used analyse the brane - world scenario in the cosmologi- 
cal context. 



A. Basic equations of the brane -world 

In Randall - Sundrum brane - world type scenarios mat- 
ter fields are confined in a three - brane embedded in a 
five - dimensional spacetime (bulk) . It is assumed that 

(5) 

the metric of this spacetime, g AB , obeys the Einstein 
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equations with a negative cosmological constant A( 5 ) 



G 



(5) 



MS ^ A (5)9ab + K (5) S (X) [~A + T Ai 



(5) 

where G)^ is the Einstein tensor, Kt^ denotes the five- 
dimensional gravitational coupling constant and Tab 
represents the energy - momentum tensor of the matter 
with the Dirac delta function reflecting the fact that mat- 
ter is confined to the spacelike hypersurface x A = \ = 
(the brane) with induced metric gAB and tension A. 

Using the Gauss - Codacci equations, the Israel junc- 
tion conditions and the Zi symmetry with respect to the 
brane the effective Einstein equations on the brane are 



(1) 



Gab — —^9ab + ^T a b + ^)S a b — £ab ■ 



(2) 



where G a b is the Einstein tensor of the induced metric 
g a b- The four - dimensional gravitational constant k and 
the cosmological constant A can be expressed in terms 
of the fundamental constants in the bulk (k(5), A(5)) and 
the brane tension A pl| . 

As mentioned in the introduction, there are two cor- 
rections to the general - relativistic equations. Firstly S a b 
represent corrections quadratic in the matter variables 
due to the form of the Gauss - Codacci equations: 

Sab — j^TTab — jT a c Tb c + ^g a b [3T cd T ca - — T 2 ] . (3) 

Secondly £ a b, corresponds to the "electric" part of the 

(5) 

five - dimensional Weyl tensor C ABCD with respect to the 
normals, ua (n A riA = 1), to the hypersurface x — 0, that 
is 



£ab = C 



(5) 
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representing the non - local effects from the free gravita- 
tional field in the bulk. The modified Einstein equations 
(^) together with the conservation of energy - momentum 
equations V a T a b = lead to a constraint on S a b and £ a b'. 



V a {£ab-K%)S ab )=0. 



(5) 



We can decompose £ a b into its ineducable parts relative 
to any timelike observers with 4 -velocity u a (u a u a — 
-I): 



£ab — 

where 



^(5) 

K 



[(u a Ub + \h a b)U + 2U( a Q b) + Tab] , 



(0) 



Q a U a = , V {ab) = Tab , V a a = , V ab U b = . (7) 

Here hi has the same form as the energy - momentum ten- 
sor of a radiation perfect fluid and for this reason is re- 
ferred to as the "dark" energy density of the Weyl fluid. 
Q a is a spatial and V a b is a spatial, symmetric and trace - 
free tensor. Q a and V a b are analogous to the usual energy 
flux vector q a and anisotropic stress tensor n a b in General 



Relativity. The constraint equation (||) leads to evolution 
equations for U and Q a , but not for Tab (see |T^|). 

The above equations correspond to the general situa- 
tion. In what follows we restrict our analysis to the case 
where £ a b = so the brane is conformally flat. Such 
bulk spacetimes admit FLRW brane - world models with 
vanishing non - local energy density (U — 0) . 



B. Scalar field dynamics on the brane 

Relative to a normal congruence of curves with tangent 
vector 



V a <5 



u a u a = -1 



(8) 



the energy - momentum tensor T^ v for a scalar field takes 
the form of a perfect fluid (See page 17 in [fl8) for details): 



with 



and 



P=\tf + V{ct>) 



V=\^ 2 ~ V{4>) . 



(9) 



(10) 



(11) 



where (j> is the momentum density of the scalar field and 
V(4>) is its potential energy. If the scalar field is not 
minimally coupled this simple representation is no longer 
valid, but it is still possible to have an imperfect fluid 
form for the energy - momentum tensor 

Substituting for p and p from ([k]) and ([ll]) into the 
energy conservation equation 

p + Q(p + p)=0 : (12) 

leads to the 1+3 form of the Klein - Gordon equation 

4> + e<p + v'(<f)) = o , (13) 

an exact ordinary differential equation for <j> once the 
potential has been specified. It is convenient to relate p 
and p by the index 7 defined by 

p=( 7 -l)p O 7=^ = -• (14) 

P P 

This index would be constant in the case of a simple one - 
component fluid, but in general will vary with time in the 
case of a scalar field: 



V 



7 = 67(7 - 2) - 2 7 ^ 



(15) 



Notice that this equation is well-defined even for (f> — > 0, 

k 
p ' 
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The dynamics of FLRW models imposed by the modi- 
fied Einstein equations are governed by the Raychaudhuri 
and Friedmann equations 



H 



H = 



H 2 



1 2 
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37 - 2 A 



(16) 
(17) 



together with the Klein Gordon equation dl3] ) above. 3 i? 
denotes the scalar curvature of the hypersurfaces orthog- 
onal to the the fluid 4 - velocity (j§) and can be expressed 
in terms of the scale factor via R = where as usual 

a 2 

k = 0,±1 determines whether the model is flat, open or 
closed. 



III. DYNAMICAL SYSTEMS ANALYSIS FOR 
EXPONENTIAL POTENTIALS 

In the following analysis we extend recent work pre- 
sented by Campus and Sopuerta to the case where the 
matter is described by a dynamical scalar field <j) with a 
self- interacting potential V(<fi) = exp(&</>), where b < 
and consider only models which have negligible non- local 
energy density U and cosmological constant A. Although 
the full dynamics are described by the equations pre- 
sented in the previous section, it is useful to re -write 
these equations in terms of a set of dimensionless expan- 
sion normalized variables. In order to get compact state 
spaces it is convenient to consider two different cases (i) 
3 i? < (k = or k = -1) and (ii) 3 R > (k = or 
k = 1). In case (i) the appropriate variables are 
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3H 2 ' 



Sli 



3 R 
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(18) 



leading via equation (|lj) to very simple expression of the 
Friedmann constraint: 

Q p + Q k + Q x = 1 ■ (19) 

We introduce a dimensionless time variable by 

(20) 

where \H\ is the absolute value of H . It follows that 

H' = -e(l+q)H, (21) 

where e is the sign of H and q is the usual deceleration 
parameter defined by 



1 d 

\H\dt 



1 a 

"iPa 



(22) 



It is clear that e — 1 corresponds to models which are 
expanding while e = — 1 correspond to models which are 
contracting. 



Using the above definitions, the dynamics for open and 
flat models are described by the following system of equa- 
tions 

Sl' k = e[(3 1 -2)(l-Sl k )+3 1 Sl x ]Sl k , (23) 
Sl' x = e[3 7 (tt A - Sl k - 1) + 2fl k }fl x . 

In case (ii) it is necessary to normalise using 



D = JH 2 + ±3i? 



(24) 



instead of the Hubble function H , so the dimensionless 
variables for this case are given by 



Q = §, si p 
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up 



3D 2 

and the Friedmann constraint is 



Sly 



K 2 p 2 

6XD 2 



Sln+Sl 



A 



1 



(25) 



(26) 



The appropriate dimensionless time derivative is defined 
via 



, _ 1 d 
= D db 



(27) 



so the dynamical equations for closed models become 



Q' = [l- 27(1 + ^a)](1-Q 2 ), 
Sl'x = 3 7 Q(^a - 1)0* ■ 



(28) 



Notice that in (|23|) , (|2q ) we have already included the 
constraint (|l9|), (|26D in order to keep the dimensionality 
of the state space as low as possible. The evolution of Sl p , 
fi„ can easily be recovered from the Friedmann equation. 



We have to add a third equation to (|23|), (28) respec- 
tively in order to describe the dynamics of the scalar field 
(j>. It turns out, that the equation of state parameter 7 
instead of <fi, is a the preferred coordinate, since it is both 
bounded by causality requirements (0 < 7 < 2) and is 
dimensionless. 

For open models (case (i)), we find using ( |l5| ) that the 
evolution of 7 is determined by 

i = e\/37(7 - 2) h/F + e sgn(4>)by/l - Sl k - Sh] (29) 

where sgn(cb) is the sign of </>. We can see from ( p9f ) that 
we only need to consider the case > 0, since the case 
ct> < can be recovered from the former by time reversal; 
simultaneously changing the sign of <f> and H results in 
an overall change of sign of 7'. From (|23]), we see that 
this transformation also changes the sign of and £1^, 
which means that <j) —> —ch corresponds to time reversal 

T — > — T. 

In the closed case, we find that 



7' = \/3t"(7 - 2)[\/3tQ + sgn(j>)bJ 1 - SI 



(30) 
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and again we can easily see from (^) and that (j> — > 
— ^ corresponds to time reversal r — * — r. Therefore in 
the following analysis we restrict ourselves to the case 
<f> > o. 

So in summary, the resulting dynamical equations we 
have to analyse are 

7 , = e v /3^( 7 -2)[V3^ + e6v / l-^-^A] , (31) 
n' fc - e[(3 7 - 2)(1 - n fc ) + 3 7 fl A ]Jl fc , 
n' A =e[3 7 (r! A -ft fe -i) + 2r! fc ]f! A 

for the open case and 

i = V37(7 ~ 2) [\/3tQ + - "a] , (32) 

Q' = [i-^7(i + ^a)](i-Q 2 ), 

h' x = 3 7 Q(ft A - 1)0 A 

for the closed case. 

Notice that in the open sector, the cosmological con- 
stant-like subset 7 = 0, the stiff- matter -subset 7 = 2, 
the flat subset Qk = 0, the GR subset fl\ — and the 
vacuum subset fi^ + ^a = 1 are invariant sets. Analo- 
gously, the closed sector has the invariant sets 7 = 0, 2, 
the flat subset Q = ±1, the GR subset fl\ = and the 
vacuum subset il\ = 1. 

IV. ANALYSIS OF THE DYNAMICAL SYSTEM 

In order to analyse the dynamical systems ([5l]) and 
(|32|), we will use for the most part the standard method 
of linearising the dynamical equations around any equi- 
librium points. One can easily see that because of the 
7' -equation, the linearisation matrix (the Jacobian), is 
not well-defined for certain values of (7, Ofe, Q\) if b ^ 0. 
Indeed, if b 7^ and 7 = or b ^ and fife + ^a = 1 
= 1), the Jacobian diverges. This means that we 
cannot linearise the dynamical equations in the neigh- 
bourhood of these points. Instead, we will have to con- 
sider the full non - linear equations and study the be- 
haviour of small perturbations away from these problem- 
atic equilibrium points. 

Furthermore the fact that we are dealing with a non- 
linear system is also important in cases, where the dy- 
namical equations can be linearised. Whenever the linear 
terms are not dominant, the peculiar behaviour of non- 
linear systems becomes visible. We then have to analyse 
the dynamical equations with additional caution. If for 
example an eigenvalue of the linearised system vanishes 
at an equilibrium point, which means that the first or- 
der terms of the dynamical equation vanishes, we have to 
study the higher order terms in a perturbative analysis 
around that point. 

Notice that the dynamical systems ([n]) and (|3^) match 
at the fc = plane (the flat subspace). If we analyse 
critical points that lie in the flat subspace by studying 



small perturbations out of that surface, i.e. where ( |3l| ) 
differs from (|32]), it is necessary to do the analysis in 
both coordinate systems and check, whether the results 
agree. If they don't, this means that small perturbations 
around the critical point will evolve differently, depend- 
ing on whether they enter the closed or the open sector. 
In our analysis we have checked this carefully. In only 

one case (F+ ) the higher order terms were dominant 
and of different sign, depending on which sector was en- 
tered. This will be commented on below. In all the other 
cases, the behaviour of the perturbations did not depend 
on the sector they entered |2^] . 



A. Models with Non- positive Spatial Curvature 

The dynamical system (pi]) has 5 hyperbolic critical 

b 2 /3 

points corresponding to the flat FLRW universe F + ' 
with 7 = f and a(t) = t 2 ^ 2 §|; thc flat FLRW uni- 
verse F 2 with stiff matter and a(t) = t 1 / 3 ; the Milne 
universe M 2 with stiff matter and a(t) = t; a flat non- 
general- relativistic model m 2 with 7 = 2, which has been 

discussed in 0; and a set of universe models X^! 3 (b) with 
7 = 2/3 and curvature fi& = 1 — 2/b 2 depending on the 
value of b. The critical points, their coordinates in state 
space and their eigenvalues is given in Table I below. In 
Table II we give the corresponding eigenvectors. Notice 

that F+ ^ only occurs for < b 2 < 6, and X+ 3 (b) for 
b 2 > 2 only. Both only occur in the expanding sector 
e = 1. 

We further find the non - hyperbolic critical points M® 
and m®(Cl\), where the former describes the Milne uni- 
verse with 7 = and thc latter is a line of non -general - 
relativistic critical points with 7 = including the flat 
FLRW model with constant energy density. The Ja- 
cobian of the dynamical system ( |3l| ) at both M e ° and 
m®(fl\) diverges for 5^0, which means that the dy- 
namical system cannot be linearised around these points 
for arbitrary positive values of b 2 . 

For b — 0, the Jacobian is well-defined for both 
M e ° and m®(Q\), and the eigenvalues and correspond- 
ing eigenvectors are given in Tables I, II. For fl\ ^ I, the 
Jacobian around m®(£l\) can still be evaluated even for 
b 0. The results from taking the limits 7, ilk —> are 
included in Table I. 

If Cl\ = 1 and b 2 > 0, we have to look at the 
non - linear system ( |3l] ) and study small perturbations 
x(r),j/(r), z(t) about the equilibrium point m^(fl\ = I). 
That means we evaluate ( |3l|) at (7, flk, £l\) = (x, y, I — z), 
where < x(t),u(t), z(t) <§; 1. 

Up to first order, ( |3l| ) becomes 

x' = —6ex — 2\IH>b\[x\J z — y) , 
V = -2ey , 
z' = -2ey . 
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TABLE I: This table gives the coordinates and eigenvalues of the critical points with non- positive spatial curvature. We have 
defined tjj = v/jj — 3. 

Model Coordinates Eigenvalues 

m^(fi A / 1) for b = (0,0,n A ) -2e(3,l,0) 

m°(n A ^ 1) /or (0,0, fix) (oo,-2e,0) 

F 2 (2,0,0) (6e + V&b, 4e, -6e) 

Me (2,1,0) 2e(3,-2,-5) 

m 2 (2,0,1) 26(3,5,3) 

F+ 2/S (£,0,0) (£-3,b 2 -2,-6 2 ) 

X?/ 3 (6) (1,1-^,0) (-1-V, -1 + ^,-2) 



For y — z or b — 0, this can be solved to give 

a; = xoe~ 6cT , 
2/ = Vae 
z = z + y a e~ 2eT . 

For y ^ z and 6^0, we find 

x = 3b 2 z a (T - t ) 2 , 
y = y a e~ 2eT , 
z = z Q + y e~ 2eT , 

where xo, yo & n( A are positive constants of integration 
(initial values). 

This demonstrates, that for 6^0, m5_(OA = 1) is 
an unstable saddle point (to be precise, it is a line of 
saddles, since z is stationary), whereas the contracting 
model m°_(fl\ = 1) remains a source (actually a line of 
sources) for all b. 

In the same way, we analyse the nature of the 
point M e °. Here, we find that the perturbed system 
(7, ri/c, fix) = (x, 1 — y, z) becomes up to first order 

x = — 6ex — 2\fib\fx\J y — z , 

V' = 2q/ , 
z = 2ez . 

If z = y or b = 0, this can be solved to give 

x = x e~ 6eT , 

V = Voe 
z = z e 

where of course yo = z a in the case y = z. 

For z ^ y and b 7^ 0, the perturbations behave like 

x = 3b 2 (y — z )(ee eT + const) 2 , 
z = z e 2eT . 

Notice that the Friedmann constraint ( |l9| ) translates 
into < y — z < 1, therefore in particular y — z > 0. 

Since we are looking at small perturbations < x(t) <C 
1, we conclude that the constant of integration is positive 



for e = — 1, whereas const < for e = 1. This shows, that 
M e ° is a saddle for all values of b. It should be realised 
that for b 7^ 0, the dynamics around the expanding model 
arc reflecting the non-linearity of the system: x is 
decreasing until e T = —const, but then the system is 
repelled since x increases for —const < e T < 00. 

Table III below summarizes the dynamical character of 
the critical points with non - positive spatial curvature. 

Note that if one of the eigenvalues in Table I vanishes, 
this just means that the lowest order terms vanish. This 
is very different to what happens in the case of systems of 
linear differential equations, where vanishing eigenvalues 
indicate lines of critical points. For a non -linear system 
to contain a line of critical points, we need small per- 
turbations away from the critical point to be stationary 
in one direction. We found an example of that case in 
the analysis of m®(fl\) above. Since we are dealing with 
a non -linear system here, vanishing eigenvalues only in- 
dicate that it is not sufficient to look at the linearised 
equations, i.e. to study the eigenvalues and eigenvectors 
of the Jacobian. Instead, we carried out a perturbative 
analysis of the kind that we have presented above when 
analyzing M e ° and m®(fl\). Including the higher order 
contributions, we obtained the results presented in Ta- 
bles III and VI. 



B. Models with Positive Spatial Curvature 

We now analyse the dynamical system ( p2| ) , which de- 
scribes flat or closed models. In addition to the criti- 
cal points corresponding to the flat models F+ F 2 , 
m®(Cl\) and m 2 , we find a maximally non -general - 
relativistic Einstein universe - like model E 1 / 3 with 7 = 
1/3, k = 1 and H = 0, which for 6 = degenerates into a 
line (hyperbola) of critical points E(£l\) with 76 [5, §]• 
For < b 2 < 2, we also find the model X 1 , (b), which 
has 7 = 2/3 and Q = ^fb/y/2. A summary of all these 
critical points, their coordinates in state space and their 
eigenvalues can be found in Table IV below. The corre- 
sponding eigenvectors are given in Table V. 

In order to analyse the nature of m®(fl\), we have to 
study the first order perturbations around the point, as 
shown in detail in the previous section. We find here, that 
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TABLE II: Eigenvalues and eigenvectors of the critical points with non-positive spatial curvature 3 R < 0. We have defined 
ip = \/ yz ~ 3- Notice that the first eigenvector of second row is (— 6e — V3b^J~^^- , 0, 3efl\ — 1)) — ► (1,0,0) as 7 — > for 

Qa ¥= i- 



Model 



m°(fi A 



Eigenvalues 



Eigenvectors 



/ 1) /or 6 = 
/ 1) /or 6 / 
F 2 

X 2/3 (6) 



e(-6,-2,0) 
(00, -2e,0) 
(6e + %/66, 4e, -6e) 
2e(3,5,3) 
2e(3,-2,-5) 



3,& 2 



2,-6 2 



(2,o,n A (i-n A )),(o,i,-n A ),(o,o,i) 
(i,o,o),(o,-i,n A ), (0,0,1) 

(1,0,0), (0,1,0), (0,0,1) 
(1,0,0), (0,1,-1), (0,0,1) a 
(1,0,0), (0,1,0), (0,1,-1) 

(1,0,0), -2), 4 + 1,0), (6 2 (1 - f ),0,3( 



1)) 



(-1-^,-1 + ^,-2) 



-2& 2 



</>,0), (-|6 2 ,1 + V,0), (-2,3(1 - A), -3(1 - A)) 



a We can use any two linearly independent vectors in the fij, = 0- 
plane as first and third eigenvectors; we have chosen the ones above 
for convenience. 



the behaviour of the perturbed system ag rees t o leading 



order with the results obtained in section IV A 



Notice that the Jacobian of the dynamical system 
around the equilibrium point corresponding to the static 
model E 1 / 3 is not well-defined for b =/= however we 
can extract the eigenvalues using a limiting procedure. 
The results are given in Table IV. Since the third eigen- 
value vanishes, it is not sufficient to analyse the linearised 
equations at that point, unless we remain in the invari- 
ant subset fl\ = 1, i n which E 1 ^ 3 appears as a saddle 
with the two eigenvectors given in Table V. A study of 
the full perturbed equations around E 1 / 3 confirms, that 
this point is a saddle however it has a more complicated 
behaviour out of the Cl\ - plane; to second order, the sys- 
tem will be repelled in the fix - direction in the expanding 
sector (Q > 0) , while in the collapsing sector (Q < 0), 
it is attracted in that direction. 

The dynamical character of all the critical points with 
non - negative spatial curvature is summarized in Table 
VI. 



V. THE STRUCTURE OF THE STATE SPACE 

With the information we have obtained in the last 
section about the equilibrium points of the dynamical 
system, we can now analyse the structure of the state 
space. The whole state space is obtained by matching 
the dynamical systems ([H ) and (|3^) . It consists of three 
pieces corresponding to the collapsing open models, the 
closed models and the expanding open models, matched 
together at the corresponding flat boundaries. 

As mentioned above, the state space has the invariant 
subsets 7 = 0, 2; fl k = Q = ±1; flxMx = and 
fi fe + fi A) A = 1. 

The bottom planes in Figures 1-9 correspond to GR, 
while top planes represent the vacuum boundaries. 

Since the eigenvalues of the equilibrium points in gen- 
eral depend on the parameter b, the nature of these criti- 
cal points varies with b. Furthermore, some of the critical 



points {F + , X + ) will also move in state space as b 
increases, since their coordinates are also functions of b. 
We also find critical points, which only appear for b = 
and then disappear, as we increase b 2 . 

For certain values of b, some of the eigenvalues pass 
through zero. At these values of b, the state space is 
"torn" in that it undergoes topological changes. In con- 
trast to the model that fij analysed, the appearance of 
such bifurcations in our model does not coincide with the 
appearance of lines of critical points. This a consequence 
of the fact, that in our inflationary model, the critical 
points are moving in state space. Bifurcations appear, 
when two of the equilibrium points merge, which occurs 
at b 2 — 0, 2, 6. These are the same parameter values as 
in GR [|l0| pT| . This was to be expected, since the only 
points moving in state space and causing these bifurca- 
tions, are the ones corresponding to the general relativis- 

tic models ^ 3 , X^ 3 . For this reason the occurrence of 
bifurcations is restricted to the GR-subspace. Therefore 
the dyna mics o f these bifurcations will only be discussed 
in section 



VA 



Since the collapsing open sector remains unchanged 
with varying parameter value b, we do not include that 
sector in the different portraits of state space. The whole 
state space can be obtained by matching Fig. 1 to Fig. 
2 -Fig. 9. 



A. The GR-subspace 

As discussed above, the subset fl\,fl\ = is the in- 
variant submanifold of the full state space correspond- 
ing to GR. In this section, we will discuss the stability 
of the general relativistic equilibrium points within the 
GR-subspace. In the next section, we will discuss the 
brane - world - modifications of these general relativistic 
results due to the additional degrees of freedom tlx, 
and also the additional non - general relativistic equilib- 
rium points. 

It is worth mentioning that although the GR- 
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TABLE III: Dynamical character of the critical points with non - positive spatial curvature. 

Model b = <b 2 < 2 b 2 = 2 2 < & 2 < § 6 2 = | f < b 2 < 6 b 2 = 6 b 2 > 6 

7n^(?Tx) line of sinks line of saddles line of saddles line of saddles line of saddles line of saddles line of saddles line of saddles 

m ( L(tl\) line of sources line of sources line of sources line of sources line of sources line of sources line of sources line of sources 

M° saddle saddle saddle saddle saddle saddle saddle saddle 

M 2 saddle saddle saddle saddle saddle saddle saddle saddle 

F+ saddle saddle saddle saddle saddle saddle saddle saddle 

Fi. saddle saddle saddle saddle saddle saddle saddle saddle 

m+ source source source source source source source source 

raL sink sink sink sink sink sink sink sink 

b 2 /3 

F + line of sinks sink saddle a saddle saddle saddle saddle 

X 2 ^ 3 (6) — - saddle a sink sink spiral sink spiral sink spiral sink 



"Notice that this point is an attractor for all open or flat models, 
while it is a repeller for all closed models. 



submanifold has been discussed in detail by Burd and 
Barrow and Halliwell pH , their analysis is some- 
what incomplete, since the state space they considered 
was non- compact. The variables they chose to describe 
the inflationary dynamics of homogeneous isotropic mod- 
els were ft, a'. Here a relates to the scale factor a(r) by 
a(r) = e a ( T \ and the potential has been absorbed into 
the time derivative by rescaling the time variable t — ► r 
by 4- — ' = V^ 1 ^ 2 -^- We find that these coordinates 
relate to our expansion normalized variables via the fol- 
lowing transformations: 



(«') 



^2 



1 



3(2-7)(l-fifc) ' 



2-7 

for 7^2, Q/. ^ 1. We can see that 

Q 

and 



oo f2fe — > 1 ^ P — > 



oo O 7 



(33) 



(34) 



This means, that we have compactified the state space 
by mapping ft G [0, oo] — ► 7 G [0, 2] and a' G [0, 00] — > 
flk G [0, 1]. In this way we have extended the work that 
has been done on inflationary models with exponential 
potentials in the general relativistic context. We find 

2/3 fe 2 /3 

that Xjl , F, correspond to the critical points I, II in 
p0| . In addition, we obtain the new critical points F 2 , 
and M 2 corresponding to a FLRW universe with stiff 
matter and a Milne universe with 7 = or stiff matter 
7 = 2 respectively. 

Notice that we also find the new equilibrium point F e °. 
The reason is that we are using 7 to describe the dynam- 
ics of the scalar field <j), which essentially is a function 
of ft 2 instead of ft Therefore our 7 -equation is homo- 
geneous in (j), whereas the dynamical equation for cf> is 
inhomogeneous in ft 

It might seem surprising, that we find a critical point 
at = 0, whereas there is no equilibrium point at <f> = 



in flO|| . The reason for this is, that we are describing dif- 
ferent physical quantities; although the equation of state 
does not change as <f) — > (7 = 0), 4> does change as 

4>^o ,(ftft' ^0). 

The collapsing FLRW model with stiff matter found in 
this analysis is of particular interest, since for all values 
of b, F 2 is the future attractor for all collapsing open and 
some of the closed models within the GR-subspace. F® 
is the past attractor for the whole collapsing open and 
parts of the closed sector for all b. 

Notice that the dynamics of the collapsing open sector 
does not depend on the parameter value b. The dynam- 
ics of this sector is constrained by the fact, that for all 
models in this sector, the flat FLRW model with constant 
energy density is the past attractor, and the flat FLRW 
model with stiff matter is the future attractor (see FIG 

I)- 

Having said that, we will now discuss the more com- 
plicated dynamics of the closed and the expanding open 
models for the different ranges of the parameter value b. 

At b = 0, we find a bifurcation of the state space, 
since for this value of b, the equilibrium points E x l z and 

X+ 3 (b), as well as F° and ^ coincide. All expanding 
open models are attracted to F®. All closed models will 
evolve into either F 2 or F9, depending on the initial 
conditions. At 7 = 2/3, the Einstein universe appears as 
an unstable saddle point (see FIG 2). 

As b 2 increases, the Einstein universe disappears. In- 
stead, we find the saddle point X + (b), which moves in 
Q- direction as b 2 is increases, remaining a saddle until it 
reaches the flat sector Q = 1 at b 2 = 2. There it merges 



with expanding FLRW model F 

V/3 



b^/3 



For < b 2 < 2, 



Fi is a sink moving along the 7 -axis. All open and 
some of the closed models are attracted to this solution 
(see FIG 3). 



At b 2 = 2, the models AT+ /3 (6) and F" + /J merge, which 
causes a bifurcation of the state space. The nature of the 
two merging critical points is significantly changed at this 
value of b. The flat solution AT+ /3 (6 = -y/2) = f! /3 is 



b 2 /3 
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TABLE IV: Coordinates and eigenvalues of the critical points with non- negative spatial curvature. We have defined \ = 



8-36 2 



<t>= V^I + 3f2 A + I 



Model Coordinates Eigenvalues 

m e (fi A / 1) /or 6 = (0, 1, fi A ) -2e(3,l,0) 

m°(n A ^ I) /or 6 / (0,1, Q x ) (oo,-2e,0) 

Fe (2,e,0) (6e + n/66, 4e, -6e) 

m? (2,e,lj 2e(3,5,3) 

F^ /3 (£.1,0) ( ^_ 3 ,^-2,-6 2 ) 

X f( fe ) (I. -73-0) (^-X,^+X,V2b) 

E 1/3 (|,0,1) (V5,-V5,0) 

E(«a)M& = (aoToI)' '^) (t^T^-T^T^ ) 



an attractor for all open or flat models, but a repeller for 
all closed models (see FIG 4). This behaviour has been 
observed by (l(J. At this value of b 2 , the two merging 
critical points swap their nature: for all b 2 > 2, X+ 3 (b) 

b 2 /3 

will be a sink, while F + will be a saddle (see FIGS 
5-9). 

b 2 /3 

As b is further increased (2 < 6 2 < 6), F + moves 
further along the 7- axis. It is now a saddle for all mod- 
els. All flat models are attracted, whereas all open or 
closed models are repelled. X+ 3 (b) has now entered the 
open sector and moves further to the fifc = 1 boundary. 
xV 3 {b) is a node sink for all 2 < b 2 < 8/3 and a spiral 
sink for 8/3 < b 2 < 00. For all b 2 > 2, F 2 is the future 

attractor for all closed models and X^ 3 (b) the future 
attractor for all expanding open models (see FIGS 5 - 7) . 

At b 2 = 6, we find another bifurcation of the state 

b 2 /3 

space. The two points F + and F 2 merge. This turns 
F 2 from a source into a saddle (see FIG. 8). For b 2 > 6, 
all open or closed models are still repelled from F 2 , but 
all flat models are now attracted. In the limit b 2 — ► 00, 
the future attractor for the open sector, X^/ 3 (b) ap- 
proaches the vacuum solution fl p — > with 7 = 2/3 (see 
FIG 9). 



B. Higher dimensional effects 

The question we will now discuss is how the behaviour 
of the dynamical system describing GR is changed within 
the brane- world context. In particular, do the additional 
degree of freedom (l\ change the stability of the gen- 
eral relativistic models, and are there new non -general 
relativistic stable equilibrium points? 

We will first answer the second question. In addition 
to the general relativistic points, we have found the non- 
general - relativistic equilibrium points m 2 and the lines 
m®(£l\). We can solve the Friedmann equation at these 
points to determine their behaviour in detail. At mj, the 
energy density is constant p(t) = po, and we can easily 



integrate the Friedmann equation to find 

a(t) ~ * , (35) 

where A = n 2 p (l + fj) behaves like a modified cosmo- 
logical constant. 

For m 2 ^, the scale factor is given by 

a(t) = (t-t BB ) 1 / (i (t + t BB ) 1 / 6 . (36) 

where t BB = 1/V6k 2 A is the Big Bang time. As 
explained in detail in this model corresponds to 
the Binetruy- Deffayet- Langlois (BDL) solution || with 
scale factor a(t) = (t — t BB ) x / e . The solution for m?_ is 
given by time reversal and t BB should now be identified 
as the Big Crunch time. 

At b = 0, we also observe, that the general - relativistic 
Einstein universe has non - general - relativistic ana- 
logues. We find a whole line of Einstein universe - like 
static equilibrium points E(fl\) extending in the fl\ di- 
rection. For b 2 > 0, the line collapses to the non- 
general- relativistic Einstein saddle E 1 ^ 3 in the fi^ = 1 
subset, and the general relativistic model X+ 3 (b). 

We now study the dynamical character of these non- 
general- relativistic equilibrium points. We find that in 
the full state space, m 2 _ instead of F 2 is the future at- 
tractor for all collapsing open models and some of the 
closed models. This means, that within the brane -world 
context, FLRW with stiff matter changes from a stable 
solution into a unstable saddle. F^_ remains a past at- 
tractor for the collapsing open and parts of the closed 
sector, but now there is a whole line m°_(flx) of sources 
extending in £! A direction. The same applies to the cor- 
responding expanding models: the sink/saddle F9_ is a 
one -element subset of the line of sinks/saddles m5_(f2^). 
This means, that for 6 = 0, the future attractor of the 
expanding open and some of the closed models is not 
necessarily F9, but instead any of the points rr&_(Clx) 
depending on the initial conditions. For b 2 > 0, m+(S}\) 
including F° turns into a line of saddles. 

For b > the models mP^^x) represent high energy 
inflationary models with exponential potentials (which 
are too steep to inflate in GR) |l3|]. The fact that they 
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TABLE V: Eigenvalues and eigenvectors of the critical points with non-negative spatial curvature R > 0. We have defined 



8-3b 2 



X = \l - o" ' f = + 3 ^A + 1. Notice that the first eigenvector of second row is (-6e 

(1,0,0) as 7 -> for Q x + 1. 

Model Eigenvalues Eigenvectors 

m^(fi A P 1) /or 6 = e(-6, -2, 0) (2, 0, fi A (l - fi A )), (0, 1, 0), (0, 0, 1) 

m°(f7 A / 1) /or b (oo,-2e,0) (1, 0, 0), (0, 1, 0), (0, 0, 1) 

F 2 (6e + V66,4 e ,-6e) (1, 0, 0), (0, 1, 0), (0, 0, 1) 

m? 2e(3,5,3) (1,0,0), (0,1,0), (0,0,1) 

F+ 2/3 (^-3,& 2 -2,-6 2 ) ( i, ,0),(6 2 (f -2),^ + l,0),(6 2 (l-f ),0,3(f - 1)) 

X+ /3 (6) (t5^X,^+X,^26) (f,^+X,0),(§,^-x,0),(2, 5^6(1-^), 3(1-^)) 

E 1/3 (>/5,-V5,0) (VS, -3,0), (V5, 3,0),- 

E(» A ) ,/or 6 = (-j^, --^,Q) 0; » A (» A _ 1)), (jgjfo^ -n A (n A - 1)), (1,0, -1(1 + n A ) 2 ) 

a We can use any two linearly independent vectors in the Q/, = 0- 
plane as first and third eigenvectors; we have chosen the ones above 
for convenience. 



are all unstable (saddle points) reflects the fact that steep 
inflation ends naturally, since as the energy drops below 
the brane tension the condition for inflation no longer 
holds. 

Finally let us consider whether the stability of the most 

b 2 /3 2 /3 . 

interesting equilibrium points F + and is changed 
by the higher - dimensional degrees of freedom. From Ta- 
bles V, VI, we can see that the third eigenvalues (corre- 
sponding to eigenvectors pointing out of the GR - plane) 

are negative for all b =^ 0. That means, if X^ 3 was 

a sink in GR, it will remain a sink in the brane -world 
scenario. 



VI. DISCUSSION AND CONCLUSION 

In this paper we extend recent work by Campos and 
Sopuerta [j to the case where the matter is described by 
a dynamical scalar field 4> with an exponential potential. 
By using expansion normalised variables which compact- 
ifies the state space we built on earlier results due to Burd 
and Barrow jl(| and Halliwell jllj] for the case of GR and 
explored the effects induced by higher dimensions in the 
brane -world scenario. 

As in J?j] we obtain the equilibrium point corresponding 
to the BDL model (m-t) Q which dominates the dynam- 
ics at high energies (near the Big Bang and Big Crunch) , 
where the extra - dimension effects become dominant sup- 
porting the claim that this solution is a generic feature 
of the state space of more general cosmological models in 
the brane -world scenario. 

We emphasise again that here, unlike in the analysis 
by Campos and Sopuerta M, 7 is a dynamical variable. 



Fixing 7 to be a constant corresponds to looking at the 
7 = const, slices of the full state space. This obviously 
only makes sense for the invariant sets 7 = and 7 = 2. 
The important point is, that even if we want to anal- 
yse the dynamical character of the de Sitter and Milne 
models in the 7 = 0- plane, we have to bear in mind the 
dynamical character of 7. Unlike in 0, we have to study 
perturbations away from the plane. This makes these 
models much more interesting in the presence of an ex- 
ponential potential, since the dynamics are not reduced 
to the plane 7 = 0. In fact, unless b = 0, we find that the 
plane 7 = is unstable. Small perturbations out of that 
plane will in general be enhanced, i.e. even for initial 
conditions with negligible 7, the system will in general 
evolve towards 7 = 2/3 or 7 = 2. Notice that orbits con- 
fined to the 7 = plane evolve towards the expanding de 
Sitter models m5j_(f2 A ) or the contracting Milne universe 
M" for all values of b. 

Finally we note that we did not find any new bifur- 
cations in this simple brane - world scenario because we 
consider only the case hi = 0. In the next paper in this se- 
ries |Q we will analyse both the effects of the non- local 
energy density Li on the FLRW brane -world dynamics 
and look at homogeneous and anisotropic models with 
an exponential potential. 
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FIG. 1: State space for the collapsing open FLRW models, e = — 1, 3 R < 0. The bottom plane Q\ = corresponds to general 
relativity. The top surface fifc + Q.\ — 1 corresponds to vacuum fi p = 0. The equilibrium points M°, Af£, F!?, ml, ml(f2,\) 
describe the Milne universe with 7 = 0/ stiff matter, the fiat FLRW model with stiff matter, the non-general-relativistic BDL 
model with stiff matter and a line of non-general-relativistic models with constant energy density (including the flat FLRW 
model with 7 = 0). We are only drawing the trajectories in the invariant planes, from which the whole dynamics can be 
deduced. The structure of this part of the whole state space does not change with the parameter value b. The full state space 
can be obtained by matching this section to FIG. 2 - FIG. 9. 




FIG. 2: State space for the FLRW models with non-negative spatial curvature 3 R > 0, e = ±1 (on the left) and the expanding 
FLRW models with non-positive spatial curvature, 3 R < 0, e = 1 (on the right) for b = (a bifurcation). In the left part of 
the figure, which describes the closed models, the plane Q = differentiates between the expanding sector Q > 0, e = 1 and 
the collapsing sector Q < 0, e — — 1. As in FIG. 1, the bottom plane corresponds to GR, whereas the top surfaces represent 
the vacuum solutions. We only give the trajectories on the invariant planes, from which the whole dynamics can be deduced. 
The critical points m^, m®(Q\), E(Q\) correspond to the flat FLRW model with stiff matter, the non-general-relativistic 
BDL model with stiff matter, a line of non-general relativistic model with constant energy density (including the flat FLRW 
model with 7 = 0) and a line of static Einstein universes with 1/3 < 7 < 2/3. 
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FIG. 3: State space for the FLRW models with < b 2 < 2. The line of Einstein universes in the closed sector of the state space 

1/3 2/3 

has collapsed into the non-general-relativistic Einstein model E ' and the general-relativistic model X + , which is expanding 

and moving towards the expanding flat subspace Q = 1, (e = 1) as b 2 is increasing. The equilibrium point ^ 3 corresponds 
to a flat FLRW model with 7 = b 2 /3. See the captions of FIG. 1, FIG. 2 for more details. 




FIG. 4: State space for the FLRW models with b 2 = 2 (a bifurcation). The equilibrium point X^ 3 has reached the flat 

fe 2 /3 

subspace, where it merges with F + /J . This causes the bifurcation; the nature of the two critical points will be swapped as they 

are moving on (see FIG. 5-9). All open and flat models are attracted to this point X + — F + , whereas all closed models 
are repelled. See the captions of FIG. 1, FIG. 2 for more details. 
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FIG. 5: State space for the FLRW models with 2 < b 2 < 8/3. X 2 J 3 has entered the open sector and turned into a node sink. 

f> 2 /3 

F, is now an unstable saddle for all models. See the captions of FIG. 1, FIG. 2 for more details. 




FIG. 6: State space for the FLRW models with b 2 = 8/3. The node sink X^/ 3 is turning into a spiral sink; the structure of 
the state space has not changed. See the captions of FIG. 1, FIG. 2 for more details. 
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FIG. 8: State space for the FLRW models with b 2 = 6 (a bifurcation). The equilibrium point ^ 3 has merged with the 
equilibrium point F+, which causes the bifurcation. The topology of the state space is changing at this value of b, since the 
point F+ is turning from an attractor in 7-direction into a repeller in that direction (compare to FIG. 7, FIG. 9). See the 
caption of FIG. 1, FIG. 2 for more details. 



